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The X-ray spectra of luminous Seyfert 1 galaxies often appear to be reflection dominated. In a number of Narrow Line 
Seyfert 1 (NLS1) galaxies and galactic black holes in the very high state, the variability of the continuum and of the iron 
line are decoupled, the reflected component being often much less variable than the continuum. These properties have been 
interpreted as effects of gravitational light bending. In this framework, we present detailed Monte-Carlo simulations of 
the reflection continuum in the Kerr metric. These calculations confirrm that the spectra and variability behaviour of these 
sources can be reproduced by the light bending model. As an alternative to the light bending model, we show that similar 
observational properties are expected from radiation pressure dominated discs subject to violent clumping instabilities and, 
as a result, have a highly inhomogeneous two-phase structure. In this model, most of the observed spectral and variability 
features originate from the complex geometrical structure of the inner regions of near-Eddington accretion flows and are 
therefore a signature of accretion physics rather than general relativity. 
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1 Introduction 



The spectra of luminous Seyfert galaxies are very well de- 
scribed by photoionized and strongly relativistically blurred 
reflection models (Fabian et al. 2004, 2005; Crummy et al. 
2006; Porquet 2006). In these sources the primary contin- 
uum often appears to be strongly suppressed. In a number of 
Narrow Line Seyfert 1 (NLS1) galaxies and galactic black 
holes in the very high state, the variabilities of the contin- 
uum and of the iron line are decoupled, in apparent contra- 
diction with the predictions of simple disc reflection mod- 
els (see, e.g. Miniutti, Fabian & Miller 2004; Fabian et al. 
2004). In particular, the monitoring of Seyfert galaxies in- 
dicates that the reflection flux can be weakly variable when 
the primary emission changes dramatically (Papadakis et al. 
2002). Moreover in at least two AGN, MGC 6-30-15 (Mini- 
utti et al 2003) and NGC4051 (Ponti et al. 2006) and one 
X-ray binary (XTE J 1650-500, Rossi et al. 2005) the reflec- 
tion flux is correlated to the primary emission at low fluxes 
and saturates at higher fluxes. In this paper we discuss two 
possible models for these puzzling properties. In the light 
bending model of Miniutti & Fabian (2004) they are inter- 
preted in terms of general relativistic effects, while, in the 
clumpy disc model of Merloni et al. (2006), these properties 
would be the result of the complex structure of near Edding- 
ton accretion flows. 
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2 Light bending model 



In this model, the active coronal region(s) illuminating the 
disc are idealised as a ring source at some height above, 
and corotating with the accretion disc. When the source is 
close enough to the black hole, the primary component is 
strongly suppressed leading to reflection dominated spectra. 
Moreover, as shown by Miniutti & Fabian (2004), fluctua- 
tions in the height of the source can lead to strong variabil- 
ity in the primary component with little variability in the 
reflected flux, as observed. These results were recently con- 
firmed by Suebsuwong et al. (2006) which present calcula- 
tions of the light bending model spectra. These calculations 
improve upon the previous works by fully computing broad 
band angle dependent reflection spectra and primary emis- 
sion as a function of ring radius, p, and height h, including 
the effects of multiple reflection. Fig.^show some sample 
spectra obtained when varying the ring radius while keep- 
ing a constant height above the disc (h = 2 R g ) . These 
spectra show that the smaller the radius, the stronger the re- 
flection component. Indeed when the primary X-ray source 
is located at less than a few R g from the black hole horizon, 
the light bending effects are very strong and tend to beam 
the primary radiation toward the disc, leading to strongly re- 
flection dominated spectra. In these extreme cases the irra- 
diation is concentrated in the central part of the disc. Then 
the same light bending effects make it difficult for the re- 
flected photons to escape. A significant fraction of them re- 
turns to the disc where they can be reflected again and so 
on. The reflection spectrum is then made up of the sum of 
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the multiply reflected spectra, which make it significantly 
different above 10 keV from what is expected from single 
reflection models (as shown in Fig.Q These effects should 
be taken into account when fitting the spectra of extreme 
sources such as MGC 6-30-15 with broad band instruments 
such as BeppoSAX, INTEGRAL or Suzaku. 

Figure|2]shows the dependence of the observed reflected 
and primary luminosity upon the height h and radius p of 
the ring source. Fluctuations of the geometrical parameters 
of the source lead to changes in the reflected and primary 
flux that are not always correlated. Indeed, as can be seen 
on this figure, depending on the parameter regime, the re- 
flected and primary components can also be anti-correlated 
or nearly independent. When the source height changes at 
constant radius and as long as p s < 5R g , its track in this 
plane can be described according to three regimes: i) at low 
fluxes (or low source height) the reflected and primary flux 
are correlated, ii) at higher fluxes the reflection saturates at 
an almost constant value while the primary can change by a 
factor larger than 2, iii) at even higher fluxes the reflection 
component is weakly anti-correlated with the primary emis- 
sion. This behaviour is described in more details by Miniutti 
and Fabian (2004). As shown by these authors many prop- 
erties of the variability of Seyfert galaxies and black hole 
binaries can be understood in terms of fluctuations of the 
source height. In particular the curious behaviour of MGC 
6-30-15, NGC4051, XTE J1650-500, where the reflection 
flux is correlated to the primary emission at low fluxes and 
saturates at higher fluxes, is in qualitative agreement with 
the predictions of this model. Fig|2]enables us to investigate 
further the model parameter space. If the radius is larger 
than ~ 5i? s the variability induced by change in the height 
is much too weak (< 2) to account for the variability ob- 
served in most accreting black holes. Let now consider the 
effects of changes in the source radius at constant height. 
At small source heights (h < 5R g ), the overall trend is that 
the reflection and primary emission are weakly correlated: 
the reflected flux changes by at most 50 % when the pri- 
mary flux increases by more than one order of magnitude 
which might be in qualitative agreement with some obser- 
vations but is inconsistent with the strong non-linear corre- 
lation observed for instance in the low state of NGC405 1 . 
At higher source heights, the reflected and primary fluxes 
become anti-correlated, which is not observed. The slope 
of the anti-correlation increases with h. At h ~ 10R g we 
could observe large variations of the reflection component 
at constant primary flux. These results suggest that if the 
light bending model is to be the correct interpretation of the 
observations, the driver of the variability should be h while 
the source ring radius has to be nearly constant and small 

(< 5Rg). 

3 Inhomogeneous accretion flows 

Black Hole systems, either in X-ray binaries or in active 
galactic nuclei, when shining at luminosities close to the Ed- 
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Fig. 1 Light bending model spectra (light blue lines) for 
a ring height of 2 R g and various radii p = 2, 5, 10 and 20 
R g , from bottom to top respectively. The thick black lines 
show the spectra obtained when multiple reflections, due to 
radiation returning to the disc, are neglected. The effects 
of multiple reflections are very significant above 10 keV . 
The dotted curves show the Newtonian model (PEXRAV) 
providing the best fit of the ligth bending spectra in the 2-30 
keV band. The reflected spectra are for neutral material and 
standard abundances. The primary emission consists in an e- 
folded power law with photon index T = 2 and E c = 200 
keV. The inclination angle is 30°. The spectra are reflection 
dominated when the ring radius is low, as can be seen from 
the reflection coefficients R > 1 shown on the figure. 
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Fig. 2 The total reflected energy flux as a function of the 
total direct energy flux, both are in the energy range 1-30 
keV, for an observation inclination 30° and different source 
height and radius. The thin lines represent the fixed source 
radius and the thick lines indicate the same source height. 
These curves are spline interpolations between the results 
of the light bending model spectra simulations lying at their 
intersections. 
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Fig. 3 Inhomogenous disc model: angle averaged spec- 
tra for tb ranging from 1 to 8 as indicated. The best fit 
parameter R and Y obtained when these spectra are fitted 
with PEXRAV in the 2-30 keV range are shown as well. 
The other fixed model parameters are the vertical Thomson 
optical of the disc tt = 1, the ionisation parameter of the of 
the cold clumps £ = 300, and the size of the regions where 
the plasma is heated h — 0.1H (see Merloni et al. 2006 for 
details). 
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Fig. 4 Inhomogeneous disc: the relative intensity of the 
reflection component (R, solid line), of the reprocessed UV 
(Luv, dot-dashed line) and reflected (Lr, dashed line) lu- 
minosities are plotted as functions of the emergent X-ray 
(Comptonised power-law) luminosity above 1 keV (£x) for 
a varying cloud optical depth t-q. All luminosities are renor- 
malized to the total heating rate Lh (see Merloni et al. 2006 
for details). The vertical dotted lines mark the X-ray lumi- 
nosities corresponding to tb=0.5,1,3,5 



dington limit are thought to be powered by accretion 
through geometrically thin, optically thick discs. Nonethe- 
less, according to the standard accretion disc solutions, such 
highly luminous discs should be radiation pressure domi- 
nated and therefore unstable to perturbations of both mass 
flow (Lightman & Eardley 1974) and heating rates (Shakura 
& Sunyaev 1973). Thus, it is not clear yet to what extent 
this standard solutions represent a realistic description of 
the observed systems. In recent years, both analytic works 
(Blaes and Socrates 2001, 2003) and simulations (Turner, 
Stone and Sano 2002, Turner et al. 2003, Turner 2004) have 
shown that magnetized, radiation pressure dominated ac- 
cretion discs may be in fact subject to violent clumping 
instabilities if magnetic pressure exceeds gas pressure and 
photons are able to diffuse from compressed regions. Large 
density variations may also be caused by photon bubble in- 
stabilities (Gammie 1998), which may develop into a se- 
ries of shocks propagating through the plasma (Begelman 
2001,2002; Turner at al. 2005). These instabilities may in 
turn have profound effects not only on the nature of the 
cooling mechanism of luminous discs, but also on their ob- 
servational appearance (see e.g. Ballantyne et al. 2004). 

Recently Merloni et al. (2006, hereafter M06) explored 
in detail the consequences of the hypothesis that accretion 
flows close to the Eddington rate are indeed inhomogeneous 
both in its density and heating structure. In this model we 
assume that the disc consist in a inhomogeneously heated 
plasma pervaded by small cold dense clumps. M06 present a 
detailed study of their expected X-ray spectral and variabil- 
ity, taking into account the radiative coupling between the 
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two phases. The spectrum emerging from this inhomoge- 
neous, two-phase, clumpy flow is calculated by coupling the 
non-linear Monte Carlo code of Malzac & Jourdain (2000), 
accounting for the Compton equilibrium between the cold 
and the hot phase, with the X-ray ionized reflection code 
of Ross & Fabian (1993, 2005), that accurately computes 
the reprocessed radiation in the cold phase. A crucial fac- 
tor in determining the broad band spectral properties of an 
inhomogeneous flow is the amount of cold clouds pervad- 
ing the hot plasma. In the limit of small scale optically thick 
clumps, this is quantified by the cloud optical depth tb = 
nSH, where H is the disc scale height , n is the cold clouds 
number density and S their average geometrical cross sec- 
tion. Figure|5]shows a sequence of spectra obtained by vary- 
ing tb- As the cloud covering fraction (i.e. tb) increases, the 
primary spectrum becomes softer, because of the enhanced 
cooling in the hot phase, and at the same time the reflec- 
tion/reprocessing features become more and more promi- 
nent. Reflection dominated spectra are achieved for tb > 1. 
As these spectra are supposed to be formed in the central 
part of the accretion flow, one naturally expect some rela- 
tivistic blurring that is not taken into account in the spectra 
of Fig. 

M06 also derived analytical formulae to estimate the 
luminosity of the different spectral components of the ra- 
diation emerging from the inhomogeneous accretion flow. 
Fig. |4] shows the reflected luminosity, the soft reprocessed 
luminosity and the reflection fraction as a function of the X- 
ray (Comptonized) luminosity. The fact that the reflected lu- 
minosity has a maximum, implies that large variations of the 
emergent X-ray luminosity, Lx, associated with changes in 
the cold clump integrated optical depth correspond to only 
modest variations of the reflection component, at least as 
long as Lx/ih *5 0.1. On the other hand, for low values of 
the Comptonized X-ray luminosity, the reflected luminos- 
ity correlates with Lx, while at high Lx, the trend is the 
opposite. This global behaviour of Lr vs. Lx is strikingly 
similar to that of the light bending model (see Fig|2j and re- 
produces the variability properties of the continuum and of 
the iron line in a number of Narrow Line Seyfert 1 (NLS1) 
galaxies and galactic black holes in the very high state . 

4 Light bending or inhomogenous accretion ? 

These similarities between the light bending model and in- 
homogeneous disc models lends itself to a simple geomet- 
rical explanation. In the more general framework of two- 
phase models for the X-ray spectra of accreting black holes, 
the main spectral and variability properties are determined 
by the geometry (and on the topology) of the two phases 
and, in particular, on the sky covering fraction of the cold 
phase as seen by the hot, Comptonising medium. Reflection 
dominated spectra are expected when the cold phase inter- 
cepts most of the photons coming from the hot phase. This, 
in the light bending model, is achieved via general relativis- 
tic effects, while in the inhomogeneous disc model it is a 



result of the clumpy and inhomogeneous nature of the inner 
disc. 

In principle, the relativistic blurring induced by the dif- 
ferential rotation of the inner disc should always be taken 
into account when fitting observed spectra. In the original 
light-bending model, where the illuminating source is a 
point-like source above a standard geometrically thin disc, 
the ratio of the reflected component to the power-law con- 
tinuum is determined by the same effects that determine the 
shape of the relativistic lines, while if the disc is truly inho- 
mogeneous, the two effects can be decoupled. Therefore, si- 
multaneous spectroscopic studies of relativistically blurred 
emission lines and of the broad band continuum and vari- 
ability could be effectively used to discriminate between a 
pure light bending model and a clumpy disc. Detailed pre- 
dictions for the latter, however, require the combination of 
sophisticated MHD and radiative transfer simulations. 
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